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Popular Summary 

on De^n.ber 1 8 , 1999. .he (MODIS). 

instruments including the Moderate R^solu | data^including global snow-cover 

Many geophysical products are thus 

products. Much of the western part of t ^nvpr is imoortant. MODIS snow-cover 

knowledge of the areal extent and ept o ,|,e currently-available 

products represent potential and 500-nr 

^”o m re"l"re \he global product are useful as input to 

climate models. 



Significant Findings 


Addittonally, there is a discussion of the validation activities so users can determine 
usefulness of the products, and associated errors. 
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On December IS. 1999, the Tena atell.te '“"f ^”* 5 ’ pMucB are derived 

Ure Moderate Resolution Imaging Spectroradiometer (MOmS> ^ g»I*)«“ 

from MODIS data including global snow-cover P^ducte^ .mSIDQ Distributed Artive Archive Center 

avarlablethroughtheNationa,^^^^^^^^ 

(DAAC) since September 13, 2000. MUUl^ snow ^ h^ause the MODIS products are 

Lancement of the currently-avaUable Also the snow- 

global and 500-m resolution, and have the capa 1 1 y P . . . ^ ^ generated for long-term 

:rrrrarwis'istn”sr:fn^^^ 

Tpa^S'l^SboS^ 

of snow cover from 40- 100 /o. Measure ° 0 02 - 0 10% Near-term enhancements include daily 

study Rom M.ch 6 2000. mv^ing MODIS data and 
field and aircraft measurements, is presented to show some early validation work. 

KEY WORDS: snow cover, snow mapping, remote sensing 

INTRODUCTION 

Srr„w.cove,i,rapsof«rcNodhcmHemisphereh^^^^ 

maps can be input to climate models. 

Or. December 18, 1999, the Earth 

complement of five instruments, one of which is the aLmated algorithms at 

(MODIS). MODIS data are now being used to produce snow-cover proOucts irom 
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Goddard Space Flight Center in Greenbelt, Maryland. The products are transferred to the National Snow 
Sd Ice Date Cente^ (NSmC) in Boulder, Colorado, where they are archived and distributed via the EOS 

Data Gateway (EDG). 

The MODIS snow-cover maps represent a potential improvement relative to hemispheric-s^le snow maps 
J^!aTavLwe toTy maiidy because of L improved spatial resolution and snow/cloud 
capabilities of MODIS, and the frequent global coverage. 

established nor has the accuracy of existing operational maps. The difiicu Y ® tn 

Tf^nv of I’^maps t it is not known which map is the “truth” (if any) and the techniques used to 

mapll maps are different, resulting in different 

else study of early validation efforts from a field and aircraft experiment near Keene, New Hampsh , 
March 2000, and field measurements in December 2000. 

BACKGROUND 
Instrument descriptions 

Mrais is an imaging spectroradiometer that employs a cross-track scan mirror, collecting optics, and a set 
S vid JdSr elSs to provide irnage^^ of the Earth’s surface and clouds in 36 discrete narrow 
speS"“ to 14.0 pm (Bames et .1., 1998). Key land-surface o^ec^ to 

study global vegetation and land cover, global land-surface change, ^Qg) The ’ 

See temperature and snow and ice cover on a daily or at 

spatial resolution of the MODIS instrument varies with spectral band, and ranges from 250 

nadir. 

^^^^;f~^^^S£^^signed to acquire calibrated radiances. The spectral coverage and 
radiometric res^inse of an existing multichannel instrument were 

spectral bands of the MODIS for measuring scientific parameters of ctoud and teir« g 

(Kino ill 1996^ The MAS with 50 spectral bands in the wavelength range from 0.55 to 14.2 pm, s 
Eal td a^AsI^^^^^^^ aircraft at an altitude of about 20 km. Data from MAS channe s 1- 0, 

fn the S e nea7-inLed and short-wave-infrared parts of the spectrum, are discussed in this paper. The 
M^ view^ side of nadir with an Earth swath width of 37.25 km. The 15-cm aperture spat. a 
tostotaneous f.eld-of-view is 2,5 mrad, or 50-m spatial resolution at nadir from the nominal aircraft he.g 


Thp T .andsat Enhanced Th e matic Mapper Plus (ETM ±) 

The ETM-H was launched on April 1 5, 1 999, on the Landsat-7 satellite 

(http://landsat,gsfc.nasa.gov/project/satellite,html). The ETM-. h^ 8 discrete 

12 5 urn and the spatial resolution ranges from 15 m in the panchromatic band, to 60 m in the thernia 
,„fa,Ttend All rf <he other bands have 30-m resolution, ETM+ data a«essed »bjor»e products 
and ordered from the USGS EROS Data Center in Sioux Falls, SD from the following Web addr ^ 
hftp.;^;ln.cr.u^ ETM. data provide a h*8h-r^«l“t.on ver at 

can be used to compare with the MODIS and operational snow-cover products. But ETM+ data are o y 
acquired once ever^ 1 6 days and are therefore not frequent enough for mapping changing snow-cover 

conditions operationally. 


Snow Maps 

NQAA Operational snow maps 
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National Environmental Satellite, Data, and Information Service (NESDIS). The Satellite Analysis Branch 
of NOAA’s NESDIS began to generate Northern Hemisphere Weekly Snow and Ice Cover analysis charts 
derived from NOAA’s GOES and POES visible satellite imager in November 1966. Maps were manually 
constructed and the spatial resolution of the charts was 190 km. However, since 1997, a new Interactive 
Multi-Sensor Snow and Ice Mapping System (IMS) is used by analysts to produce products daily at a 
spatial resolution of about 25 km, and utilizes a variety of satellite data to produce the maps (Ramsay, 
1998). NOAA also produces a daily product, developed by automated techniques, that uses visible, near- 
infrared and passive-microwave data to map snow cover, and agrees in 85% of the cases studied, with the 
IMS product (Romanov ct al., 2000). 

The National Operational Hydrologic Remote Sensing Center (NOHRSC) maps. NOHRSC snow-cover 
maps, generated by National Weather Service NOHRSC hydrologists are distributed electronically in near 
real time, to local, state and federal users during the snow season (Carroll, 1995). The NOHRSC 1-km 
maps are generated primarily from the NOAA polar-orbiting satellites and the Geostationary Orbiting 
Environmental Satellite (GOES) satellites to develop daily digital maps depicting the areal extent of snow 
cover for the coterminous United States, and Alaska, and portions of southern Canada. 


Description of the MODES snow-mapping p roducts and algorithrn s 

The MODIS snow products are provided as a sequence of products beginning with a swath product, and 
progressing, through spatial and temporal transformations, to an 8-day global-gridded product. For a more 
complete d^cription, see Riggs et al. (2000). The MODIS maps provide global, daily coverage at 500-m 
and, at ~5.6-km resolution at the equator. The climate-modeling grid (CMG) products are at l/20°th 
resolution which is ~5.6-km resolution. Quality-assessment (QA) information is included in the products. 
In the future, percent snow cover or fractional snow cover in each pixel will be provided along with daily 
snow albedo in the 500-m products. Fractional snow cover for each pixel (from 40-100%) is already 
available in the CMG products. Because cloudcover often precludes the acquisition of snow cover from 
visible and near-infrared sensors, the daily maps are composited, and 8-day composite products are 
available. The climate modeling grid (CMG) product is specifically designed for use by climate modelers. 
It provides a global view of the Earth’s snow cover in a geographical projection (Figure 1). 

The automated MODIS snow-mapping algorithm uses at-satellite reflectances in MODIS bands 4 (0 545- 
0.565 pm) and 6 (1.628-1.652 pm) to calculate the normalized difference snow index (NDSI) (Hall et al, 
1995): 


NDSI = hand 4 - band 6 

band 4 + band 6 [1] 

A pixel in a non-densely-forested region will be mapped as snow if the NDSI is >0.4 and reflectance in 
MODIS band 2 (0.841-0.876 pm) is >11%. However, if the MODIS band 4 reflectance is <10%, then the 
pixel will not be mapped as snow even if the other criteria are met. This prevents pixels containing very 
dark targets such as black spruce forests, from being mapped as snow. This is required because very low 
reflectances cause the denominator in the NDSI to be quite small, and only small increases in the visible 
wavelengths are required to make the NDSI value high enough to classify a pixel, erroneously, as snow. 

Changes that occur in the spectra of a forest stand as it becomes snow covered can be exploited to map 
snow cover in forests. The primary change in reflectance occurs in the visible wavelengths as snow has a 
much higher visible reflectance than soil, leaves or bark. A fundamental change that snow cover causes in 
the spectral response of a forest, which can be used in a global algorithm, is that the reflectance in the 
visible will often increase with respect to the near-infrared reflectance. This behavior is captured in the 
normalized difference vegetation index (NDVI), as snow will tend to lower the NDVI. MODIS bands I 
(0.620-0.670 pm) and 2 (0.841-0.876 pm) are used to calculate the NDVI. The NDVI and NDSI are used 
together to improve snow mapping in dense forests. If the NDVI “0.1, the pixel may be mapped as snow 
even if the NDSI is <0.4 (Klein et aL, 1998). 
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The .n., eduction cfa“.h=n.al—;ne.g, see 

Oh October 3, 200 1, hae .ntproved the snow 

found in nutny parts of the globe ,n the W^ron To tntprove 

cover were confus.orr with cloudcover, ' f we instituted a thermal mask in the 

the snow-mapping results, and remove most of the spunous .r.,:- ’ moDIS infrared bands 31 (10.78- 

fall of 2001 which has greatly reduced the errors of commission. Us g estimate 

11 77-12 27 uih) a split-window technique (Key et al., IW ), 

11.28 pm) and ( • ’ r * i "o the nlxel will not he mapped as snow. 

dramatically. 

rrSnt::r:3reist=w.r« 

specific to each snow algorithm and data product, og. suspect d^ extreme viewi g 

Additional details concerning the QA may be found in Riggs et al. (2000). 

TU.unni^rlnud mask product The current version of the MODIS snow algorithm reads the 

in snow-covered areas. 

The snow algorithm is not applied when the surface viewed is in darknras. 

comes from the cloud mask product, and is defined where the solar zenith angle _85 

Great Lakes). 

m the product. Inputs to the swath snow-cover product are. the MODIS (I^vel I B) fiance a 

'==S=SS3S 

mask product. The constraints are applied in the order listed. 

Dally and 8-day compome lile prodacs. The secoud praduct in the sequCTce » 
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™ny obs«v.,io„s ao<,u,r=d during a da, (Figure ^ S » ^ ^'y »0- 

con.pos.te n»xin.un. snow^ver product .s pro<l“^ ^ “'J’ '* of U.e daily CMC 
m resolution products, 8-day per.ods beg.n from the f“' ^ ^^Srnding^^^ CMC con.,»site 

early products and are no longer produced.) 

In selecting the compositing technique for the 8-day i^fon on the daily tile 

use mmrimum snow cover. In other words, 'f “77*“ cover .s a more useful 

?rremHi^:mrora7.a“:«^^^^ 

b^^|!L\rwr„Tnr=^^ 

minimizes cloud cover. 

ftr/fy and «-<4.y con,porr re cImate-mMng grid product. The eight-day composite climate-modeling gnd 

(CMG) products are available to download from ^ ^ and from the National Snow 

through NSIDC, and is presented m a roduct and binning the 500-m cell 

(approximately 320) to include all land ar^s, o ® ^ .j binning technique in use creates a 

observations into the ~5.6-km spatial resolution of the CMG c^s a fractional cloud data 

fractional snow cover array based on the number of snow observations in a gna ceii a ir^ 

indication of the quality of the 500-m rtata w,th.n j p„,s), using a special 

can be difficult because ma '^n« ™a mSum sl;w e^tent’map by defining. 

only from 40-100% snow cover. 

obscuration during the period. All 8 days ot o a cell during the period This results in 

snow extent is reported as the highest fraction of snow ^ ent clouds are 

on three separate 8-day composite snow maps ranged from 0.02-0.10/.- 
Data archiving and distribution. The NSIDC is one of eight NASA 

(DAACs) and is part of the Earth Observation System Data and Information System (EOS ). 1 n 
EOSDIS utilizes the EOSDIS Core System (ECS) for data rn^ jmem « ,„d ord;r 
which facilitates online Web-based user access to data (Scharfen f ,bout the 

ami limiialiani. The cloud mask (Ackermm. « U., 1*981 as “ ^ 
snow and ice algorithms, tends to overestimate cloud cover as seen m ^ ^ ^ ^ 

Twalh produc. .he Eai. Coas. of Ihc Uniled S<a.^. ^ “ Sdtask’lJnTiguTsB 

rn™-cov\mdTmr.mThowr^^ lie f'™*! '» ■"<“ '* 
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When .he Coud mash .s used. Compaq o, .his “ 

a mom realis.ie depicion of .he ac.ua. snow L .o i.s use of MODIS 

sirtrshrLS^ 

™p"Srhrh.;?sZ Sve,,. w..h ».= e«ep.K,n of high clouds .ha. eon.ain ice 

which may often be mapped erroneously as snow. 

mapping are under investigation. 

Contemn in .he iden.,f.ca.ion of cloud ove, snow has also been ote^wed in highmievalion regions, e g , 
the Sierra Nevada in California and the Southern Alps of New Zealand. 

Anolher lim,.a..on of d,e MODIS “g^t^lC— nts " 

devdop an automated algorithm that is effective globally to map fractional snow co 

Zrv^ml er S d^^^^ 

km There have been many studies employing both passive microwave and f^hfi^sMble 

cover, and following the launch of Aqua, a great leap of progress in this regard should possib 

There are mslances where snow cover ,s depos,.ed and mete "“e SU my ™w 

after .he snow has meted, .he MODIS 8-day " tlte^ .0 be 

during that 8-day period. This is an error of omission and has not been estimatea yei, du pp 

rather minor 


CASE STUDY 
Northeastern United States 

primary sites were studied. Bretwood Golf Course and Tenant S wamp, toth 

6). Skll conditions were almost completely clear for the entire day over Keene on March 6, 

Snow conditions at Bretwood Golf Course were patchy while at Tenant ^ On 

continuous, with 3-10 cm snow depths and snow densities ranging from 350-420 g cm at both sites u 
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Figure 6 Bretwood Golf Course can be seen clearly on both the MAS image and the MAS-denved snow 
map because more snow is visible than m the surrounding areas, including Tenant Swamp, which are 
forSted and therefore the snow cover underneath is partially obscured. A spherical densiometer was used 
to measure the tree-canopy density at Tenant Swamp, Results showed that the percent of open canopy 
ranged from 37-51%. Because the tree canopy is obstructing the view of the ground from above, e 
MODIS and MAS snow maps do not show 1 00% snow cover on and near Tenant Swamp. 

A snow map was derived from MAS data using a prototype MODIS snow-mapping f ^ 

shown in Figure 6. Patchy snow cover is evident at the site of the field measurements with the Bretwood 
Golf Course showing as snow covered even though field measurements reported patchy snow Jhe 

smaller patches can L estimated as 4 x 8 m in size. The total snow-covered area at the golf course site was 
-60%. Patches range from very small: less than 1 m^ to large continuous patches 

and almost totally continuous toward the NW part of the golf course near toe river, MODIS from toe 
same day also show patchy snow near Keene. Though the snow cover at Tenant Swamp was continuous, 
according to ground measurements, because of the forests, not all of the snow on the ground is visible. 

The MODIS snow map and the NOHRSC operational snow map were compared for an 

part of Maine, all of New Hampshire and Vermont, Massachusetts, and ^rt of New York, Rhode H a 

and Connecticut on March 6, 2000 (Figure 7). The Keene, New Hampsh^, study arra is 

map Though there was general agreement between the MODIS and NOHRSC maps dunng the day of the 

field work on March 6, there is only -73% agreement in the area of New York and parts of New England 

when toe MODIS and NOHRSC maps are digitally registered. The 2 order RMS error when the 

NOHRSC image was registered to the MODIS image was 0.65 pixel. 

The primary area of disagreement between the snow maps is in ^uth-central 

NOHRSC snow map shows complete snow cover in toe Catskill Mountains, and toe MO p 

patchy snow cover. 

In south-central New York, 26 meteorological stations that are located in the area in question (see Fipire 
7) were found (NOAA, 2000). Of these stations, only three reported measurable snow on toe ground, six 
reported a trace of snow and the rest of the stations showed no snow was on the ground. The stations that 
had snow on March 6*^ tended to be located at the higher elevations; all of the stations (e^ept for o"e) that 
reported any snow were over 350 m in elevation. In this forested area, it is hard to detect all of the snow on 
the ground, especially if toe sensor is not ‘looking’ straight down (nadir). 

Figure 8 shows part of the NOAA NOHRSC and IMS operational maps for toe area. The IMS ™p^o 
sh^s snow cover in Canada, but that was deleted for Figure 8 to facilitate companson vnth the NOHRSC 
map which does not show snow in this part of Canada. Though at different resolutions, the IMS map 
shows less snow in the Catskills than does the NOHRSC map. 

ETM+ images from March 8, 2000, show some patchy snow cover Trough cloudcover. Combined 
analysis of images in visible bands of toe Sea-viewing Wide Field-of-view Sensor (S^WiFS) and MODIS, 
the conclusion is that the NOHRSC map overestimates toe snow cover, m the Catskills, on March 6, 2000, 
but the exact amount of snow cover is unknown. 


To summarize, four different data sets were studied to reveal the extent of snow cover in south-central New 
York. Unfortunately an ETM+ scene from March 8, 2000, was cloudy. A SeaWiFS image was available, 
and some snow was present as determined from analysis of the visible bands. A rather large 
cover was mapped on the NOHRSC snow maps contrasting with a rather small area map^d by the MODIS 
map for March 6, 2000. Snow was also visible on toe IMS product as seen in Fi^re 8. Since none ot the 
snow maps agree, the accuracy of the MODIS snow map is impossible to determine. 


On December 23, 2000, field work was again performed at the same sites near Keene, and, althoug t ere 
was no aircraft overflight, there was a Landsat-7 ETM+ acquisition on that day. According to groun 
measurements, the snow cover was continuous and snow depths varied from ~6-8 cm at all of the site 
The MODIS snow map was registered to a snow map that was developed from ETM+ data (using tne 
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MODIS snow-mapping algorithm as modified for use with ETM+ data). While the 30-m resolution ETM+ 
derived snow map showed patchy snow cover, the MODIS map, digitally registered to the ETM+ map, 
showed nearly continuous snow cover (Figure 9). The binary MODIS snow-mapping algorithm will map 
snow cover if ~50% of the 500-m pixel is snow covered. In this case where the snow is patchy within each 
500-m pixel, the MODIS snow map will overestimate the snow cover. 


DISCUSSION AND CONCLUSION 

A sequence of MODIS snow-cover products is presented. The swath products are mapped to the 
integerized sinusoidal grid to create the daily tile product. Eight days of the daily tile products are used to 
produce the 8-day composite tile product. These products are at 500-m resolution. The climate-modeling 
grid (CMC) product will be produced at ~5.6-km resolution and consists of daily and 8-day composite 
products. Examples of the products are shown, focusing on the site of a field and aircraft experiment from 
March 6, 2000, and field work on December 23, 2000. The MODIS snow map shows patchy snow cover 
on March 6th, as confirmed by the field measurements. However nearly complete snow cover is mapped 
near Keene on December 23rd, which is an overestimation of snow cover as compared to the ETM+- 
derived map for the same day, as confirmed by field measurements. 

Other work shows that the MODIS snow-cover maps compare favorably with current operational maps 
(Hall et al., in press), and perform better than do passive-microwave-derived snow-cover maps during the 
daytime and in the fall months in the Northern Hemisphere when wet snow causes the current passive- 
microwave snow-mapping algorithms not to map snow that is known to be there. 

Future enhancements to the MODIS snow maps include daily snow albedo (Klein et al., 2000), which 
should be available in early 2002 and, fractional snow cover at 500-m spatial resolution (Barton et al., in 
press; Kaufman et al., submitted). Future re-processing will allow all of the MODIS snow maps to be 
processed in a consistent manner beginning with the September 13, 2000 data. 
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FIGURES 

1 8-day composite global MODIS snow-cover map - December 1 8-25, 2000. Snow is white, ‘ha^ 

persisted fm all 8 days are pink and non-snow covered areas are light green; grey areas represent areas 

not imaged by MODIS due to darkness. . i • i j u 

2 . MODIS^image (left) and snow map (right) - November 3, 2000. Snow is ^7iZnd 

non-snow covered areas are green; water is either dark blue (oceans or bays) or light blue (inland 

3 ^^^Location of tiles on the integerized sinusoidal grid, 3B Daily tile snow wver product showing 
snow cover in the Black Hills, South Dakota. Snow is white, clouds are pink and non-snow-covered 
areas are green. Inland water is blue/green. 

4A Daily MODIS snow maps of North America - December 1 8-25, 2000; 4B 8-day compos 
MODIS Inow map of North America - December 18-25, 2000. Snow is white, 
for 8 days are pink and non-snow covered areas are light green; grey areas represent areas not imaged 

by MODIS due to darkness (SZA>85°). x/fr^nic 

MODIS snow map product of the Eastern United States - December 23, 2000 (left) 
map processed without the MODIS cloud mask of the same area (nght) - December 23, 2000. Snow 

white clouds are pink and non-snow-covered areas are green. Inland water is blue/green. 

MAS image (left) and snow map (right) showing field sites near Keene, New Hampshire (at arr ) 

MODIS^OHRSC difference map showing differences in the location and amouiit 
the eastern United States - March 6, 2000. The MODIS snow map was proc^sed without the MODIS 
cloud mask in this figure. Graph shows snow depths measured at meteorological stations in sou 
central New York - March 6, 2000. 

8 NOAA Interactive Multisensor Snow and Ice Mapping System (IMS) snow map 

National Hydrologic Remote Sensing Center (NOHRSC) snow map (right) from March 6, 2000 ‘he 
northeastern United States. Snow cover, shown on the IMS map in Canada, was removed in order to 

facilitate the comparison with the NOHRSC map. u;« 

9 Landsat-7 ETM+ (left) and MODIS (right) false-color images (top row) of Keene, New Hampshi re and 
the surrounding areas; ETM+ (left) and MODIS (right) snow maps of the same area (bottom row) - 
December 23, 2000. 
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